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Transport in graphene nanoribbons with an energy gap in the spectrum is considered in the 
presence of random charged impurity centers. At low carrier density, we predict and establish 
that the system exhibits a density inhomogeneity driven two dimensional metal-insulator transition 
that is in the percolation universality class. For very narrow graphene nanoribbons (with widths 
smaller than the disorder induced length-scale), we predict that there should be a dimensional 
crossover to the ID percolation universality class with observable signatures in the transport gap. 
In addition, there should be a crossover to the Boltzmann transport regime at high carrier densities. 
The measured conductivity exponent and the critical density are consistent with this percolation 
transition scenario. 
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One of the remarkable experimental findings of the 
past two years, creating a great deal of activity and con- 
troversy, has been the observation that the carrier density 
dependent dc conductivity of gated 2D graphene layers, 
while being approximately linear in density at high gate 
voltage, becomes a non-universal constant over a finite 
gate voltage range AVg around the charge neutral Dirac 
point. While this conductivity minimum plateau forma- 
tion around the charge neutrality point in 2D graphene 
is experimentally well established [1|, the actual values 
of cTmin ^ 2 e'^/h - 20 e'^/h and AVg - 1 y - 15 1/ are 
non-universal and depend on the sample quality [3, l3|. 
The minimum conductivity plateau in graphene has been 
theoretically explained [J, |5| to arise from the invari- 
able presence of unintentional random charged impuri- 
ties at (or near) the graphene-substrate interface which 
lead to inhomogeneous electron-hole puddle formation in 
the low gate voltage regime [J, S Ifl] ■ We note that dis- 
tortions of the graphene membrane and quenched rip- 
ples can also give rise to density inhomogeneities [7|, 
and there have been recent theories studying the effect 
of ripples on graphene conductivity [Sj . While we focus 
here on charged impurity induced inhomogeneities, many 
of our conclusions are only sensitive to the existence of 
the inhomogeneous density landscape (i.e. electron-hole 
"puddles"), and these do not distinguish between mecha- 
nisms (e.g. impurities, ripples) producing these puddles. 
Since graphene is a 2D semimetal (or more appropriately, 
a zero-gap 2D chiral semiconductor with electron-hole 
bands touching each other linearly at the charge neutral 
Dirac point), the conductivity becomes an approximate 
constant when the gate voltage induced chemical poten- 
tial is pinned in this electron-hole puddle region around 
the Dirac point. This inhomogeneous electron-hole pud- 
dle based theoretical understanding of the graphene min- 
imum conductivity plateau formation leads immediately 



to an important fundamental question: Are there situa- 
tions where this inhomogeneous puddle picture leads to 
a graphene 2D metal-insulator transition (2D MIT) as is 
known [9J, |lO[ to occur in 2D semiconductor systems? 

We show in this Letter that indeed, as a direct con- 
sequence of the inhomogeneous puddle formation in 
graphene, the system will manifest a 2D MIT, which 
is precisely in the same universality class as the cor- 
responding 2D MIT in electron [9| and hole 10] GaAs 
systems, provided that there is an energy gap separating 
the graphene electron and hole bands. The fundamental 
physics here is that of percolation - for usual 2D zero-gap 
graphene, percolation through the puddles is allowed at 
all gate voltages, occurring either through the electron 
puddles or the hole puddles (or through both ll|), since 
one or the other is always percolating. If there is a gap, 
however, then there should be a percolation-driven 2D 
MIT in graphene exactly as found [9|, [lO| in 2D GaAs 
based semiconductor structures. 

The easiest way to introduce an energy gap in 
graphene, which would then immediately lead to a 
percolation-induced transport gap (i.e. two separate 2D 
MIT transitions for electrons and holes), is to consider 
graphene nanoribbons instead of bulk 2D graphene. In 
this Letter, we predict and confirm experimentally that 
graphene nanoribbons exhibit a 2D MIT in the low car- 
rier density regime as a function of the applied gate volt- 
age and that this MIT is in the percolation universality 
class; furthermore, we predict theoretically that as rib- 
bons become very narrow, there should be a dimensional 
crossover to the ID universality, implying that the ob- 
served transport gap would tend to infinity as the ribbon 
width goes to zero (or in practice, becomes smaller than 
the typical size of the puddles), reflecting the ID percola- 
tion universality where metallic conduction is completely 
suppressed. We speculate that such a 2D-1D crossover 
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FIG. 1: Evidence of percolation driven metal insulator tran- 
sition in a graphene nanoribbon. Main panel shows graphene 
ribbon conductance as a function of gate voltage. Best fits 
at low density to Eq. [l] give for electrons A = 1.485, ric = 
26.7485 X 10" cnT'^ ,&" = 1.3 ± 0.2 and for holes A = 
1.755, ric = -18.5 x lO" cm-^(5'' = 1.3 ± 0.1. Best fit at 
high density to Eq. [2] gives riimp = 22 x 10^" cvcC^ . Inset 
shows the same data in a linear scale, where even by eye the 
transition from high density Boltzmann behavior to the low 
density percolation transport is visible. 
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FIG. 2: Percolation driven metal-insulator transition in three 
additional graphene samples. The left panel shows a naturally 
occurring graphene nanoribbon with dimensions W ~ 200 nm 
and L « 11 /xm and has a critical exponent &^ = 1.2±0.2. The 
center panel is the W = 24: nm sample reported in Ref. [Ij] 
which has a critical exponent of S^ = 1.3 ± 0.1, and the right 
panel is the W — 49 nm sample reported in Ref. [12| which 
has a critical exponent of 5'' = 1.6 ± 0.3. 



may have been observed in recent experiments l^, UM ' 
but more quantitative work and more data would be nec- 
essary to establish this prediction. 

In Fig. [1] we show our experimental data supporting a 
percolation driven metal-insulator transition on a natu- 
rally occurring graphene nanoribbon (that is expected to 
have smoother edges than comparable ribbons fabricated 
using the method of Ref. [I4I) of dimensions L « 1.5 fim 
and W ~ 50 nm. The conductance G — {W/L)a is 
shown as a function of applied gate voltage Vg oc n. 
We performed Quantum Hall measurements on the large 



area connected to the ribbons to confirm that we have 
a mono-layer of graphene and the details of the fabri- 
cation can be found in Ref. [3|. Best fits to Eq. [T] at 
low density give conductivity exponents S'^ = 1.3 ± 0.2 
and S'^ = 1.3 ± 0.1 (close to the theoretically expected 
value S — 4/3), and the fit to Eq. [2]at high density gives 
rtjnip = 22 X 10^° cm~^ which is consistent with measure- 
ments on similarly prepared bulk graphene samples [3|. 
Fig. [2] shows a similar analysis for a Vl^ = 200 nm sample 
fabricated at Maryland, and two of the Columbia samples 
(W = 24 nm and W = A9 nm) reported in Ref. ^IJi]. All 
three samples show the low-density percolation universal- 
ity class with critical exponents S = 1.2 zL 0.2, 1.3 ± 0.1, 
and 1.55 ± 0.3 respectively, which are similar to perco- 
lation exponents observed in 2D GaAs systems [Sl, Il0l |. 
The corresponding fit parameters for the three samples 
are WA/L = 0.008, 0.002 and 0.0087 respectively, and 
?ic = -65.08, 170.87 and 184.23 x lO^*' cm-^ respec- 
tively (the units oi A are [(10^° cm'^)* 25.8 kn]'^). 

For the 2D percolation universality class, at low den- 
sity we have 



A{n - n,y 



(1) 



where S « 4/3 is the 2D percolation critical exponent. 
For graphene ribbons, we expect two such percolation 
transitions, one for electrons and one for holes, separated 
by a "transport gap" defined as A^ = "/\/Tr{nl — n^), 

where nc is the critical density for electrons (holes) and 
7 = hvF is the graphene Fermi velocity. For larger carrier 
densities, where |£^f| ^ Ag, we expect a crossover __to a 



high-density Boltzmann transport regime where [J, |5| 
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just as for bulk graphene on a Si02 substrate (where 
riD is the charge neutrality point and riimp is the 2D 
surface impurity density of Coulomb scatterers). If we 
define ^ as the typical size of the electron or hole pud- 
dle, where below we calculate ^ self-consistently using 
the random phase screening approximation, then so long 
as the sample width W > ^, we would have 2D per- 
colation whereas li W < ^ one has ID percolation i.e. 
a chain of approximately L/^ p-n junctions. Changing 
nimp (which could be extracted from high density mobil- 
ity measurements) would also change ^ and the critical 
width for which this dimensional crossover is observed. 
Moreover, we predict that signatures of reduced dimen- 
sionality should be apparent in temperature dependent 
transport measurements. 

To date, most theories for transport in graphene 
nanoribbons consider a quasi ID, rather than the 2D 
limit. The experimental observation of 2D percolation 
in these ribbons casts strong doubt onto the relevance 
of the quasi ID theories for current graphene nanorib- 
bon experiments. As was already discussed in Ref. 14| . 



many features of the quasi ID geometry get washed out 
for W ^ 10 nm which is the case in most experiments 
on graphene nanoribbons. The length scale controlling 
the crossover from quasi ID to 2D behavior in this con- 
text may very well be determined by other independent 
parameters which are unknown at this stage such as the 
inelastic scattering length or the phase breaking length 
(both of which depend on temperature). Whether the 
transport properties in graphene nanoribbons should be 
considered using a 2D or a quasi ID Hamiltonian is at this 
stage an open question requiring further theoretical and 
experimental study. Our analyses involving data from 
two groups, our own and that of the Columbia group [12l |. 
clearly establish that depending on the value of W either 
the 2D or ID percolation universality class may apply, 
where we believe this transition to be controlled by the 
size of the electron and hole puddles induced by charged 
impurities. We can not rule out the possibility that fur- 
ther lowering of temperature would lead to quasi ID be- 
havior [l5| and the percolation-driven 2D MIT is only 
a crossover phenomenon. Although we focus on single- 
layer graphene, we note that a similar percolation tran- 
sition should also be seen in graphene bilayers, where 
since an electric field induced gap can be introduced into 
the spectrum without any confinement, the crossover to 
a quasi ID regime would not arise. We note that even 
for bulk graphene, a Boltzmann to percolation crossover 
could be induced with a magnetic field, where for small 
field and within the electron-hole puddle model, we ex- 
pect the p-n resistance to be very low justifying the Boltz- 
mann picture, whereas for large magnetic field, the p-n 
junction becomes very resistive jl6] inducing a percola- 
tion transition. This crossover may have been observed 
in recent experiments [17|. 

To reinforce the point that the 2D MIT in graphene 
nanoribbons is indeed a percolation transition and not a 
quantum crossover phenomenon, we calculate the perco- 
lation critical density rip using the non-linear screening 
argument of Efros [l8| with the basic idea being that 
the MIT occurs when inhomogeneous density fluctua- 
tions created by the charged impurities can no longer 
be screened by the carriers. This leads to Up ^ y/rii^/d 
where the random charged impurities of concentration 
nimp are assumed to be located at a distance d from 
the 2D graphene plane. Taking d '^ 1 nm and nimn ~ 
2 - 5 X 10" cm-2, typical values estimated 0, S li, i| 
from mobility measurements, we get n^ ~ 5 x 10^^ cm~^. 
This is in reasonable agreement with our experimental 



■' 2 X 10^°cm ^ for the same 
q <ti Up, and our experimental 



finding in Figs. [T] and [2l On the other hand, the quan- 
tum localization crossover density Uq can be estimated 
from the loffe-Regel criterion kpl ~ 1, where i is the 
mean-free path, to be Uq 
rii and d values. Thus, n. 
critical density agrees with the percolation critical den- 
sity, providing further support for a percolation driven 
insulating transition in graphene. 

Experimentally, one can measure three different gaps. 
In addition to the transport gap Ag discussed above, the 
temperature dependence of the conductivity minimum 
gives an activated gap Aact (we performed this mea- 
surement on the 200 nm ribbon, and found AactW^ ~ 
0.1 eV nm which is an order of magnitude smaller than 
theoretical estimates |19| of the confinement induced 
gap), and finally, Rcf. J12!| reported the source-drain 
bias required to induce conduction and found the gap 
to be orientation independent (contrary to the expec- 
tation of the quasi ID theory that is extremely sensi- 
tive to whether the edge is zig-zag or armchair). The 
connection between these three experimental gaps and 
the theoretical gap in the energy spectrum is beyond 
the scope of this work, where we focus here only on the 
transport gap Ag and predict that in the 2D regime, pro- 
vided that the impurity location is pretty much the same, 
|n^ — nj,'| oc ^riimp oc /i~^'^, where u is the high-density 
mobility in the Boltzmann regime [3, y, y, l5| • 

The 2D percolation picture presented above breaks 
down when the sample width becomes smaller than the 
typical disorder length scale. Using the self-consistent 
RPA method of Ref. 5] we can obtain an integral ex- 
pression for the potential correlation function {V{r)V{0)) 
which for experimentally relevant parameters can be ap- 
proximated by 
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Using rs = e'^ /kj, where k is the effective dielectric con- 
stant that depends on the choice of substrate, we find 
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where for z = Akpd and Ei{x) = J t ^e *^dt 
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This notation is chosen to be consistent with Ref. [a| 
where the rms density n* = 2rlninipCo{z = Ad^/nn*). 
For typical values of nimp ~ 20x 10^° cm"^ and d « 1 nni, 
we have ^ « 10 nm, which is consistent with the ex- 
perimentally observed critical width W* = 16 nm [12| . 
In the ID limit, there should be no percolation transi- 
tion, only an activated conduction, and the effective gap 
should diverge in the T ^ limit. From Eq. [U we can 
predict the dependence of the puddle size (and therefore 
the critical width of the dimensional crossover) on exper- 
imentally tunable parameters. For example, we predict 
that cleaner samples (i.e. with a larger high-density mo- 
bility) would have larger critical widths and by doping 
graphene with potassium 3, thereby changing only njmp, 
one could tune through this dimensional crossover in a 
sample of fixed width. Changing the substrate to a high- 
re material like Hf02 (assuming that riimp and d remain 
unchanged) could significantly increase the puddle size, 
in contrast to suspended graphene [20| where increasing 
Ts (which decreases the puddle size) is compensated by 
the lower nimp (which tends to increase the puddle size). 
In addition, we predict that suspended nanoribbon ex- 
periments will have a smaller transport gap due to the 
order of magnitude higher mobility, but the same critical 
exponent S; while for nanoribbons with potassium dop- 
ing of various strengths, increasing rjimp should lead to a 
larger transport gap with no change in S. 

In conclusion, we have argued theoretically and demon- 
strated experimentally that a disorder induced, density 
inhomogeneity driven percolation transition is observable 
in graphene nanoribbons. We anticipate a crossover to 
Boltzmann transport at high carrier density and a di- 
mensional crossover for sample widths that are smaller 
than the disorder induced puddle size. Several features 
of the experiment including the difference between trans- 
port and activation gaps, the large discrepancy between 
the value of the gap and that predicted by band structure 
theory as well as there being no orientation dependence 
of the gap and a critical width below which the dimen- 
sional crossover causes a divergence in the transport gap 
(for an infinite ID system in the T — > limit) are all 
explained naturally in this picture. This consistent theo- 
retical picture better captures the physics of nanoribbons 
than the quasi ID induced models that have dominated 
the literature to date. Our discovery of a percolation 
driven graphene 2D MIT also shows the close concep- 
tual connection between 2D graphene transport and 2D 
semiconductor transport, and establishes that density in- 
homogeneities dominate carrier transport in both classes 
of systems at low enough carrier densities. In bulk 2D 
graphene, which is a zero-gap semiconductor, this leads 
to the low-density minimum conductivity plateau, and 
in graphene nanoribbons (as well as in bilayer graphene 
with an electric field induced gap), where there is an en- 
ergy gap between the electron and hole bands, we get the 
same percolation induced 2D MIT familiar from 2D semi- 



conductor electron [9| and hole [lO| systems. For very 
narrow graphene nanoribbons, which are in the ID perco- 
lation regime, our theory predicts an insulating behavior 
with an effective infinite energy transport gap at T = 0, 
which may have been observed experimentally [12l [l3| , 
by virtue of the absence of a percolation transition in ID. 
A fundamental question of considerable significance that 
remains open in this context is the experimental absence 
of quantum localization 21| which may be observable at 
much lower temperatures than used experimentally so 
far. 

This work is partially supported by U.S. ONR. 
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